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Optimization  of  Low-Drag  Bodies  of  Revolution 
With  Turbulent  Boundary  Layers* 

by 

E.  W.  Geller,  J.  E.  Mercer,  R.  W.  Metcalfe,  and  E.  M.  Murman 
Flow  Research  Company,  Kent,  Washington  98031 


Abstract 

Numerical  optimization  methods  are  used  to  search  for  minimum-drag 
bodies  of  revolution  in  incompressible  flow  with  all-turbulent  boundary 
layers.  Available  methods  for  potential  flow,  boundary  layer,  drag,  and 
optimization  calculations  are  used.  Results  obtained  with  an  empirical 
drag  formula  give  a low-fineness-ratio  body  with  signif icantly  lower 
drag  than  for  traditional  shapes.  However,  when  the  drag  method  of 
Nakayama  and  Patel  is  used  as  a representative  state-of-the-art  capa- 
bility, results  show  that  the  drag  is  insensitive  to  perturbations  from 
the  shape  of  traditional  bodies,  and  only  a small  reduction  in  drag  is 
obtained. 
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Nomenclature 


Definition 

reference  area  for  the  drag  coefficient 
t h 

the  i parameter  defining  the  body  shape 
drag  coefficient 
skin-friction  drag  coefficient 

pressure  drag  coefficient 

1?  , drag  coefficient  based  on  the  body  volume 

hpU  V ' 

^ 00 

skin-friction  coefficient 
body  drag 

maximum  body  diameter 
^2 

t—  , boundary  layer  shape  factor 


, boundary  layer  planar  shape  factor 

61 

body  length 
body  radius 
U L 


V 


U V 


, Reynolds  number  based  on  the  body  length 
1/3 


, Reynolds  number  based  on  the  body  volume 


, Reynolds  number  based  on  the  planar  boundary- 


layer  momentum  thickness 

R/L  , non-dimensional  body  radius 


the  value  of  r at  x 


x . 

l 


velocity  component  in  the  x direction  at  the  edge  of 
the  boundary  layer 

undisturbed  free-stream  velocity 

boundary  layer  velocity  component  in  the  X direction 
body  volume 

)■  2 

tt  / (R/L)  d(x/I.)  , body  volume  based  on  unit  length 

o 
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Nomenclature  (Cont'd) 


Definition 


boundary  layer  velocity  component  in  the  y direction 

axial  coordinate  for  the  body 

non-dimensional  axial  coordinate  for  the  body 

Cartesian  coordinate  in  the  direction  of  the  body 
surface  (a  coordinate  for  the  boundary  layer  velocity 
field) 
th 

the  i non-dimensional  axial  coordinate  for  the  body 

Cartesian  coordinate  in  the  direction  normal  to  the 
surface  of  the  body  (a  coordinate  for  the  boundary 
layer  velocity  field) 

coefficient  in  equations  (16)  and  (17) 

coefficient  in  equations  (16)  and  (17) 

6 

2tt  J (1  - -jj)  p dy  , boundary  layer  displacement  area 


2tt  J (1  - -jj)  p dij  , boundary  layer  momentum  area 


boundary  layer  thickness 


j (1  - -p)  dy  , boundary  layer  planar  displacement 


thickness 


/ (1  - -jy)  -jj  dy  , boundary  layer  planar  momentum 

J thickness 

o 


slope  of  the  body  at  the  tail 

kinematic  viscosity  of  the  fluid 

distance  to  the  body  axis  (a  coordinate  for  the 
boundary  layer  velocity  field) 

fluid  density 
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1.  Introduction 

The  design  of  low-drag  bodies  is  a subject  that  hydrodynamicists 
and  aerodynamicists  have  studied  for  years.  It  is  clear  that  a reduction 
of  drag  for  a given  body  will  result  in  reduced  requirements  for  propul- 
sion systems  and  a greater  range  for  a given  fuel  supply.  In  recent 
years,  considerable  research  has  been  aimed  at  obtaining  extensive 
regions  of  laminar  flow  to  reduce  the  friction  drag  of  incompressible 
bodies  of  revolution.  Many  factors  affect  the  occurrence  of  transition 
from  laminar  to  turbulent  flow,  including  adverse  pressure  gradients, 
surface  temperature  and  roughness,  noise  and  turbulence  present  in 
the  background  media  or  generated  by  the  propulsion  system,  and  the  use 
of  suction  to  control  the  boundary  layer  velocity  profile.  In  the 
design  of  a "laminar  flow  body"  the  objective  is  to  use  the  effects  to 
achieve  as  large  a transition  Reynolds  number  as  possible  for  a practically 
shaped  body.  Unless  special  efforts  are  made  in  the  design,  manufacture, 
and  maintenance  of  the  body,  underwater  vehicles  will  experience  boundary 
layer  transition  near  the  nose,  and  almost  the  entire  surface  will  be 
covered  by  a turbulent  boundary  layer. 

The  principal  question  addressed  by  this  study  is  whether  significant 
drag  reduction  can  be  obtained  for  bodies  of  revolution  by  body  shaping 
if  it  is  assumed  that  the  boundary  layer  becomes  turbulent  very  near 
the  nose.  Previous  work  has  been  done  on  this  problem.  Gertler  (1950) 
conducted  a systematic  experimental  study  of  a family  of  streamlined 
bodies  at  the  David  Taylor  Model  Basin;  this  investigation  is  often 
referred  to  as  the  Series  58  study.  Early  work  in  airship  design 
considered  the  drag  of  bodies  of  revolution  (for  example.  Freeman  (1932a) 
and  Young  (1939));  however,  minimum-drag  shapes,  per  se,  were  not  studied. 
Parsons  (1972)  and  Parsons,  Goodson,  and  Goldschneid  (1974),  using  a 
boundary  layer  and  potential  flow  computer  program  to  calculate  the  drag, 
considered  the  application  of  numerical  optimization  methods  to  search 
for  the  shapes  of  low-drag  bodies.  Their  study  was  primarily  directed 
toward  designing  laminar  flow  bodies,  although  several  calculations  were 
made  for  all- turbulent  boundary  layers.  Using  a finite-difference 
boundary  layer  method,  Hess  and  James  (1975)  also  considered  the  possi- 
bility of  low-drag,  turbulent  boundary  layer  bodies  in  their  calculations 
of  drag.  The  results  of  all  the  above  studies  indicate  that  the  drag  of 
an  all-turbulent  boundary  layer  body  in  incompressible  flow  is  relatively 
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insensitive  to  body  shape,  provided  the  body  is  sufficiently  streamlined 
to  avoid  separation.  The  question  has  arisen,  however,  as  to  whether 
these  studies  are  conclusive  and  whether  some  more  unconventional  shape 
with  a lower  drag  might  exist. 

In  the  study  reported  here,  we  have  used  numerical  optimization 
methods  to  search  for  minimum-drag  bodies  of  revolution  for  incompressible 
flow  and  all-turbulent  boundary  layers.  Although,  our  method  is  similar 
in  spirit  to  that  used  by  Parsons  (1972),  it  is  different  in  many  respects. 
First,  our  optimization  algorithm  is  the  "method  of  feasible  directions," 
which  was  used  by  Hicks,  Murman,  and  Vanderplaats  (1974)  in  the  design 
of  low-drag  transonic  airfoils.  Second,  the  hydrodynamics  program 
includes  a potential  flow  method  and  an  integral  boundary  layer  method 
for  computing  the  necessary  parameters  to  evaluate  the  body  drag.  We 
decided  that  an  integral  boundary  layer  method  would  be  sufficiently 
accurate  and  would  require  less  computer  time  than  a finite-difference 
boundary  layer  method.  Third,  we  selected  a new  way  of  parameterizing 
the  body  shape,  which  specifies  body  ordinates  at  a discrete  number  of 
locations  as  the  "design  variables"  and  uses  cubic  splines  under  tension 
to  define  the  body  at  other  locations. 

In  the  remaining  sections  of  this  report,  the  formulation  of  the 
problem,  organization  of  the  computer  codes,  and  calculated  examples 
are  explained.  In  the  discussion  and  conclusions  section  we  observe  two 
principal  results  of  the  study.  First,  using  "state-of-the-art"  potential 
flow  and  boundary  layer  methods  and  drag  prediction  formulas,  we  did 
not  produce  body  shapes  with  significantly  lower  drag.  We  found  one 
shape  with  significantly  lower  drag,  but  the  drag  formula  used  is 
questionable.  Second,  the  "state-of-the-art"  boundary  layer  methods  are 
inaccurate  in  the  aft  body  region,  where  the  boundary  layer  is  thick 
compared  to  the  body  radius.  Boundary  layer  methods  are  inapplicable 
in  this  region  because  of  the  importance  of  radial  pressure  gradients 
and  the  interaction  of  the  viscous  and  inviscid  flow.  For  conventional 
bodies  with  high  fineness  ratios  and  thin  boundary  layers  over  most  of 
the  body,  the  inaccuracy  in  the  tail  region  has  little  effect  in  the 
prediction  of  the  total  drag.  However,  for  low-drag  body  shapes  with 
smaller  fineness  ratios  and  thicker  boundary  layers,  the  inaccuracy  in 
the  tail  region  is  significant.  This  observation  has  been  independently 
observed  in  the  recent  note  by  Patel  and  Guven  (1976).  Good  computational 
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raodels  for  this  region  are  required  before  additional  design  studies  are 
conducted . 

Finally,  we  note  that  the  basic  optimization  methods  used  in  the 
present  study  may  be  employed  in  other  design  problems  where  the 
accuracy  of  the  drag  calculation  is  not  a factor.  Such  studies  include 
maximizing  the  transition  Reynolds  number  by  optimizing  the  surface 
temperature  distribution  for  a given  heating  or  cooling  supply,  optim- 
izing suction  distribution  for  a given  pumping  capacity,  and  maximizing 
the  propulsion  system  efficiency. 


lb  “ 
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2.  Formulation  of  the  Problem 


2. 1 Statement  of  the  Problem 

The  problem  is  to  find  a body  of  revolution  of  prescribed  volume, 
which  has  minimum  drag  at  a given  speed  through  a prescribed  incompressible 
fluid.  We  assume  that  surface  quality  (in  other  words,  roughness) 
promotes  early  transition  of  the  boundary  layer  so  that  nearly  the  entire 
body  surface  has  a turbulent  boundary  layer.  Furthermore,  we  exclude 
bodies  that  have  boundary  layer  separation. 

1/3 

We  select  the  cube  root  of  the  volume  V as  the  reference 
length  scale  for  the  Reynolds  number  and  for  the  drag  coefficient, 
which  are  given  respectively  as 


U V 


1/3 


Rev 


(1) 


W2/3 


(2) 


With  this  selection,  the  Reynolds  number  is  fixed,  and  the  drag  coefficient 
is  to  be  minimized  during  the  search  for  an  optimum  body. 

To  make  the  problem  tractable,  we  limit  the  search  for  an  optimum 
body  to  a parameterized  family  of  shapes.  It  is  convenient  to  use  the 
body  length  L as  the  length  scale  for  defining  the  family  non-dimens ionally , 
and  we  write 


R/L  = f(X/L;a..  ,a~  , • • • a ) , X = 0,L  , 

1 Z n 

R/L  = 0,X=0,X  = L, 

where  a.  , i = l,n  are  free  parameters.  The  two  length  scales, 

1/3  1 

V , are  related  by  the  volume  formula 
L 1 

V = 7T  j R2  dX  = L3tt  j (R/L)2  d (X/L)  = L3V  , 
o o 

where  V is  the  body  volume  for  unit  length  given  by 


1 

V = tt  / (R/L)2  d (X/L)  . 


(3) 


L and 


. u> 


(5) 


o 
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Since  V is  fixed  and  V may  be  computed  from  the  shape  according  to 
the  preceding  formula,  there  is  only  one  independent  length  scale  for 
the  problem,  and  the  other  is  directly  related  through  the  formula 

L = (V/V)1/3  (6) 


when  the  body  shape  is  specified.  In  presenting  the  bodies,  we  plot 
all  results  using  the  variables  X/L  and  R/L  so  that  the  bodies  are 
all  shown  with  unit  length.  The  fatter  bodies,  however,  will  be 
physically  shorter  than  the  slender  bodies  for  a fixed  volume. 

The  problem  under  study  may  be  formally  stated  in  the  following 
way:  Given  a body  described  by  the  equation 

R/L  = f (X/L ;a, ,a0 , • • *a  ) , X = 0,L 
i.  L n 


R=0,X=0,X  = L, 


(3) 


where  L is  the  length  of  the  body,  R is  the  radius,  X is  the  axial 
coordinate,  and  a^  , i = l,n  are  free  parameters,  and  given  a method 
of  computing  the  drag  coefficient  of  the  body  for  a specified  Reynolds 
number  and  boundary  layer  transition  point 


CD  = F(R/L,Rev)  , 

find  the  values  of  a^  that  minimize  the  drag  coefficient,  subject  to 
various  constraints.  For  the  present  problem  the  constraints  selected 
are 

(1)  That  the  body  have  a nonnegative  radius 


R(X;a^ ,a9 , *an)  > 0 , 0 < X < L , 


(2)  That  the  Reynolds  number  constant 

A/3 


U V 


Rev  = 


(7) 


v 


constant  , 


(8) 
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(3) 


which  implies  for  a given  speed  and  a given  fluid  (in  other  words, 
a given  and  V ) that  the  volume  be  constant 

L 

V(a^  ,a2 , • • *an)  = 7T  J R^  dX  = constant  , (9) 

o 

That  the  boundary  layer  be  turbulent  (in  other  words,  that  the 
transition  point  be  at  the  nose*) , and 


(4)  That  boundary  layer  separation  does  not  occur  on  the  body 


Cf(X;a1,a9 , *an)  > 0 , X = 0,L  . 


(10) 


2. 2 Optimization  Method 

The  optimization  method  is  basically  a numerical  search  procedure 
in  which  we  start  with  an  initial  body  shape,  and,  by  systematically 
altering  the  values  of  the  design  variables  a^  , we  reduce  the  objective 
function,  the  drag  D , subject  to  the  imposed  constraints,  until  ve 
reach  a minimum  value.  This  design  method  is  a direct  method  because  the 
body  shape  is  specified,  and  the  flow  field  and  drag  are  calculated. 

Various  numerical  optimization  methods  are  available.  In  this  study, 
we  used  the  "method  of  feasible  directions,"  which  was  implemented  by 
Vanderplaats  (1973)  and  which  has  been  successfully  applied  to  problems 
involving  the  optimal  shaping  of  transonic  airfoils.  A detailed 
explanation  of  the  methodology  is  given  by  Vanderplaats,  Hicks,  and  Murman 
(1975).  The  procedure  is  briefly  outlined  here. 

With  the  initial  values  of  the  design  variables  a^  and  the  objective 
function  D specified,  we  can  calculate  the  gradient  of  D by  using 
the  following  one-sided  finite  difference  formulae 


3D  D(ai,***a1+Aai, •••an)  - D(at , • • •a1 , • • • afl) 

Sa^  Aa^ 


(ID 


Assuming  that  no  constraints  are  violated,  we  then  proceed  with  our 
search  for  a minimum  drag  shape  in  the  direction  of  steepest  descent; 


In  this  study  we  approximated  this  condition  by  fixing  transition  at 
five  percent  of  the  length  downstream  from  the  nose. 


r 

♦ 


i 
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that  is,  in  the  direction  opposite  to  the  gradient  direction.  Predetermined 
step  sizes  in  the  space  are  used,  and  three  steps  are  taken.  The 

final  step  size  is  selected  with  the  conjugate  gradient  method  of  Fletcher 
and  Reeves  (1964).  After  this  step  search  is  completed,  the  gradient  of 
D is  recalculated,  and  another  search  is  done.  Assuming  that  no  constraints 
are  violated  in  the  process,  this  search  procedure  is  continued  until 
the  change  in  drag  in  three  successive  searches  is  less  than  a specified 
convergence  criteria. 

If  one  or  more  constraints  are  violated  (or  nearly  violated),  then 
the  gradient  of  the  constraint  is  also  calculated  by  finite-difference 
formulae,  and  the  search  proceeds  in  the  direction  that  reduces  the  drag 
and  avoids  violating  the  constraint;  that  is,  the  "feasible  direction." 

Again  the  search  proceeds  as  above  (without  using  the  conjugate  gradient 
method)  until  convergence  is  achieved  or  until  it  is  determined  that  no 
feasible  solution  exists.  Clearly,  the  final  solution  in  either  case 
may  not  be  an  absolute  minimum.  It  may  be  a local  minimum  or,  if  the 
design  space  is  "flat,"  the  calculations  may  terminate  because  the  conver- 
gence criterion  is  satisfied  in  three  successive  searches.  In  either 
case,  however,  if  a significant  reduction  in  drag  is  achieved,  the  method 
may  be  considered  successful  from  an  engineering  viewpoint. 

2. 3 Body  Parameterization 

For  the  optimization  method  to  be  efficient,  the  body  must  be 
described  by  as  few  parameters  as  possible.  On  the  other  hand,  for 
the  optimization  method  to  be  general,  the  body  parameterization  method 
chosen  must  admit  as  diverse  a set  of  shapes  as  possible.  Thus,  the 
selection  of  the  body  shape  parameterization  is  an  important  part  of  the 
optimization  calculation. 

Several  possible  body  parameterization  methods  were  considered, 
but  we  actually  tried  only  one.  k'e  list  those  considered  since  they  may 
be  useful  for  future  optimization  studies. 

The  simplest  parameterization  is  a polvnominal  family 

N 

r2  = E a x"  , 0 < x < 1 (12) 

1 ° 

like  the  one  used  in  the  Series  58  study  (Gertler,  1950).  Hicks,  Murman, 
and  Vanderplaats  (1974)  used  a similar  family.  A more  versatile 


The  use  of  different  formulae  on  the  front  and  back  part  of  the  body 
allows  each  region  to  be  modified  more  independently.  However,  continuity 
of  r and  one  or  more  of  its  derivatives  is  enforced  at  the  match  point 
x^  . Parsons  (1972)  and  Vanderplaats , Hicks,  and  Murman  (1975)  used 
this  procedure.  Another  possible  method  is  to  use  super-ellipses 


where  the  functions  f (x)  are  selected  to  have  certain  desirable 

n 

properties.  Hicks  and  Vanderplaats  (1977)  used  this  method  very 
successfully  in  the  design  of  low-drag  airfoils. 

The  method  used  in  this  study  is  to  choose  the  body  radii 
at  a number  of  fixed  locations  x^  as  the  parameters  and  then  to  use 
a cubic-spline-under-tension  interpolation  method  reported  by  Cline 
(1974)  to  describe  the  body  at  intermediate  values  of  x . This  method 
has  the  advantage  that  the  body  ordinates  themselves  become  the  design 
variables.  The  spline  method  requires  that  either  the  first  or  second 
derivatives  of  r be  specified  at  the  end  point.  Thus,  for  the  tail 
point,  we  made  the  body  slope  9 a design  variable.  At  the  nose, 
however,  the  slope  is  generally  infinite.  In  this  region,  we  used  the 
formula 


(16) 
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and  matched  the  value  of  r and  its  first  two  derivatives  at  x^  with 
the  spline  interpolation  for  the  remainder  of  the  body.  In  our  early 
calculations  (reported  as  Calculation  No.  1 in  section  4)  we  used  the 
simpler  formula 

al 

r = ax  , 0 < x < x^  (17) 

and  matched  only  r and  its  first  derivative  at  x^  . The  matching 
with  the  spline  determines  a and  3 . The  exponent  a^  , however, 
remains  free  as  a design  variable,  and  a variety  of  nose  bluntnesses 
are  therefore  possible.  To  summarize,  for  the  present  study,  the  design 
variables  are  the  following: 

nose  bluntness  factor 
body  radii 
tail  angle 

In  the  calculated  examples,  N is  5. 

Two  variations  of  the  above  method  were  considered  but  not  tried. 

One  variation  is  to  allow  the  abscissa  locations  x^  to  be  design 
variables  also.  The  other  variation  is  to  consider  the  perturbation  from 
a fixed  basic  body  shape.  In  this  procedure,  the  initial  body  would  be 
the  basic  body  shape,  and  the  perturbed  radii  would  be  the  design  variables. 
If  the  basic  shape  has  a blunt  nose,  the  special  nose  shape  function, 
equation  (16),  would  not  be  required. 

2. 4 Hydrodynamic  Model 

Once  a body  shape  is  specified,  a hydrodynamic  calculation  must  be 
performed  to  compute  the  drag.  The  drag  arises  from  the  viscous  forces 
and  may  be  separated  into  two  components:  (a)  the  shear  stress  or 

tangential  force  on  the  body,  known  as  the  friction  drag;  and  (b)  the 
normal  stress  or  pressure  force,  known  as  the  pressure  drag.  The  sum 
of  the  pressure  drag  and  the  friction  drag  is  the  total  or  profile  drag. 

This  is  the  quantity  we  want  to  minimize. 

The  standard  procedure  for  calculating  the  drag  is  to  first  calculate 
the  inviscid  flow  with  a potential  flow  method  and,  from  this  result,  to 


r(xi)  , i = 1,N 

6 (18) 
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obtain  the  pressure  distribution  on  the  body.  Next,  with  a boundary 
layer  method,  the  viscous  layer  is  calculated  up  to  the  separation  point. 
Finally,  the  drag  is  evaluated  from  Young's  (1939)  or  Granville's  (1953) 
formula. 

In  the  present  study,  several  assumptions  were  made  in  assembling 
the  hydrodynamics  model.  We  assumed  that  for  low-drag  bodies,  separation 
would  not  occur  ahead  of  the  tail,  and  therefore  our  viscous  flow  method 
did  not  need  to  treat  separation.  We  also  assumed  that  since  we  are 
principally  interested  in  an  integrated  force  — the  drag,  an  integral 
boundary  layer  method  would  suffice  and  would  be  substantially  faster 
than  a finite-difference  method.  The  results  published  by  Hess  and  James 
(1975)  support  this  assumption.  The  third  assumption  was  that  influence 
of  the  viscous  flow  region  in  the  immediate  vicinity  of  the  tail  would 
not  substantially  affect  the  overall  drag,  and,  hence,  a boundary  layer 
method  could  be  used  along  the  entire  length  of  the  body  even  though  it 
is  incorrect  in  the  tail  region.  This  assumption  was  based  on  published 
methods  for  evaluating  the  drag  with  Young's  and  Granville's  formula. 
Unfortunately,  this  assumption  turned  out  to  be  the  weak  link  in  the 
hydrodynamics  model,  as  we  discuss  later  in  this  report.  For  low-drag 
bodies  with  thick  boundary  layers  at  the  aft  end,  the  breakdown  of  the 
boundary  layer  assumption  is  significant.  The  result  has  recently  been 
reported  by  Patel  and  Guven  (1976). 


In  the  next  section,  we  describe  the  actual  boundary  layer,  potential 
flow,  and  drag  formulae  that  were  used.  We  note  here,  however,  that 
two  separate  programs  were  employed.  Computer  Program  No.  1 employs  a 
potential  flow  method  that  Frank  Woodward*  developed  by  using  a source 
distribution  on  the  body  surface  according  to  the  method  of  Hess  and  Smith 
(1966).  The  boundary  layer  method  was  developed  by  Frank  Dvorak*  using 
the  Nash-Hicks  (1968)  integral  method.  Calculations  performed  with  this 
program  produced  an  interesting  shape.  The  result,  however,  is  questionable 
because  of  the  thick  boundary  layer  on  the  aft  end.  We  therefore  decided 
to  use  a program  by  Nakayama  and  Patel  (1974),  which  is  well  documented 
and  verified  for  traditional  streamline  body  shapes.  This  program, 
referenced  as  Computer  Program  No.  2,  was  obtained  directly  from 
Professor  Patel. 

* Current  Address:  Analytical  Methods,  Inc.,  Bellevue,  Washington. 


i Flow  Research  Report  No.  79 

April  1977 

-U- 

3.  The  Computer  Programs 

A system  for  computing  a body  shape  with  minimum  drag,  according 
to  the  preceding  discussion,  has  several  well-delineated  tasks: 

(1)  Parametrically  defining  - ;dy  shape, 

(2)  Calculating  the  potential  flow, 

(3)  Calculating  the  boundary  layer, 

(4)  Coupling  the  potential  flow  and  boundary  layer  calculations, 

(5)  Calculating  the  drag,  and 

(6)  Searching  for  the  parameter  values  giving  minimum  drag. 

Two  computer  programs  were  developed  by  assembling  various  codes 
available  in  the  literature  and  at  Flow  Research  Company.  The  flow 
diagrams  shown  in  figures  1 and  2 show  our  choices  for  the  tasks  listed 
above . 

3. 1 Computer  Program  No.  1 

Computer  Program  No.  1 uses  the  method  of  Hess  and  Smith  (1966)  for 
'*  the  potential  flow  calculation.  The  flow  field  is  constructed  with  sources 

distributed  on  the  body  surface.  The  governing  equation,  a Fredholm 
integral  equation  of  the  second  kind,  enforces  the  desired  normal 
velocity  at  the  surface  and  is  solved  numerically  by  approximating  the 
surface  source  distribution  with  constant-strength  planar  panels. 

In  Computer  Program  No.  1,  the  method  of  Nash  and  Hicks  (1968)  is 
used  to  calculate  the  growth  of  the  turbulent  boundary  layer.  It  is 
based  on  the  momentum  and  the  moment  of  momentum  integral  relations 
for  the  boundary  layer  and  on  Cole's  (1956)  velocity  profile  family. 

The  potential  flow  and  boundary  layer  calculations  are  coupled 
according  to  the  traditional  weak  interaction  approach.  A potential 
flow  is  calculated  to  obtain  a pressure  distribution,  which  is  then 
used  for  input  to  the  boundary  layer  growth  calculation.  The  potential- 
flow  calculation  is  repeated  with  the  boundary  layer  growth  effects 
modeled  by  ejecting  fluid  from  the  body  surface.  The  ejection  velocity 
(in  other  words,  the  boundary  condition  for  the  potential  flow  calculation) 
is  calculated  according  to  the  relation  developed  in  appendix  A.  The 
new  pressure  distribution  is  used  to  repeat  the  boundary  layer  calculation. 
If  greater  accuracy  is  desired,  the  cycle  may  be  repeated;  however,  for 
our  purposes,  additional  iteration  was  not  necessary. 
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To  check  the  coupled  boundary  layer  and  potential  flow  calculation, 
we  analyzed  a body  of  revolution  for  which  extensive  experimental  data 
were  obtained  by  Yi  and  Przirembel  (1973).  The  body  that  was  tested  was 
a long  cylinder  followed  by  an  elliptical  aftbody.  This  body  could 
not  be  directly  modeled  in  our  calibration.  Instead,  we  used  a shortened 
body  with  three  segments  of  equal  length:  the  first  segment  was  a 4 to  1 

elliptical  forebody,  the  second  segment  was  a cylinder,  and  the  third 
segment  was  a 4 to  1 elliptical  aftbody.  The  initial  development  of 
the  boundary  layer  was  impossible  to  model  since  the  experimental  body 
extended  upstream  of  the  tunnel  test  section.  Therefore,  both  the 
displacement  and  momentum  thicknesses  over  the  elliptical  aftbody  have 
smaller  initial  values  in  the  computed  results  than  the  experimental 
data  indicates.  The  shape  factor,  however,  is  almost  correct,  so  this 
parameter,  along  with  pressure  coefficient,  provided  an  adequate  means 
of  evaluation.  Figure  3 shows  the  pressure  coefficient  and  shape  factor 
comparisons.  For  this  calculation,  separation  was  predicted,  and  the 
extrapolation  of  the  normal  ejection  velocity  beyond  that  point  may 
explain  the  deviation  of  the  calculated  pressure  from  that  measured  near 
the  tail.  The  shape  factor  agrees  well  in  the  region  of  decelerating 
flow. 

To  calculate  the  drag  in  Computer  Program  No.  1,  we  first  compute 
the  friction  drag  by  integrating  the  surface  shear  stress  from  the  boundary 
layer  calculation,  and  we  then  use  an  empirical  relation  attributable  to 
Young  (1939).  Young  noted  that  theoretically  derived  pressure  drags  were 
approximately  proportional  to  the  thickness  ratio  and  to  the  total  drag, 
and  he  proposed  the  following  formula: 


CD  = 0.4(D/L)Cd  , (19) 

P 

where  is  the  total  drag  coefficient,  is  the  pressure  drag 

coefficient,  0 is  the  maximum  body  diameter^  and  L is  the  body 
length.  Using  this  relation  we  compute  the  total  drag  from  the  calculated 
friction  drag  according  to 


C 


D 


1 - 0. 4( D/L)  ’ 


(20) 


i 
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where  is  the  skin  friction  drag  coefficient.  The  advantage  of 

using  this  simple  formula  is  that  the  friction  drag  is  accurately 
calculated  by  our  method  while  the  pressure  drag  is  difficult  to  obtain 
by  direct  integration. 

To  assess  the  accuracy  of  the  drag  calculation  method  in  Computer 
Program  No.  1,  we  used  it  to  compute  the  drag  for  six  of  the  Series  58 
bodies.'*'  The  comparison  with  the  measured  drag  given  in  table  1 and 
figure  4 shows  that  the  procedure  is  adequate  for  our  purposes  when  the 
fineness  ratio  is  seven  or  above.  For  the  boundary  layer  calculations 
no  iterations  were  made  with  the  potential  flow  calculation  to  account 
for  displacement  thickness.  To  assure  that  the  inaccurate  drag  cal- 
culations at  low  fineness  ratio  were  not  the  result  of  this  neglect,  we 
repeated  the  calculation  for  the  Model  4154  (fineness  ratio  4) , using 
one  iteration.  There  was  little  change  in  the  result. 


Table  1 

L/D  Cn  , Drag  Coefficient  at  a 

v 6 

Series  58  Finess  Reynolds  No.  Re^  = 2 x 10 

Model  No.  Ratio  

Measured*- 
Gertler  (1950) 


4154 

4 

0.0197 

0.0165 

4155 

5 

0.0196 

0.0176 

4156 

6 

0.0197 

0.0186 

4165 

7 

0.0206 

0.0202 

4158 

8 

0.0221 

0.0216 

4159 

10 

0.0198 

0.0195 

Computed  by  eq.  (1) 


*~We  refer  frequently  to  the  Series  58  bodies,  which  are  a family  of  shapes 
tested  at  the  David  Taylor  Model  Basin  and  reported  by  Gertler  (1950). 

2 

A trip  strip  five  percent  of  the  length  downstream  from  the  nose  provided 
a turbulent  boundary  layer. 


Error  Factor  Based  on  Table 
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The  standard  procedure  for  calculating  drag  based  on  the  well-known 
Squire-Young  formula  was  discarded  after  making  seme  comparisons  with 
experiment  as  described  in  the  following.  The  formula  given  by  Young 
(1939)  for  this  procedure  is 


where  A is  the  reference  area  for  the  drag  coefficient  CD  , and  where 
the  following  are  evaluated  at  the  tail  of  the  body: 

A 2 = the  boundary  layer  momentum  area, 

U = the  velocity  at  the  edge  of  the  boundary  layer,  and 
H = the  boundary  layer  shape  factor  . 

Using  the  pressure  distribution  measured  in  the  Langley  20-foot  propeller 
research  tunnel  for  the  Airship  Akron  and  reported  by  Freeman  (1932b)  , we 
calculated  the  boundary  layer  growth  and  then  computed  the  drag,  first 
according  to  equation  (20)  and  then  according  to  equation  (21).  The 
comparison  with  the  measured  drag  reported  by  Freeman  (1932b)  is  given 
in  table  2.  For  a second  comparison,  we  analyzed  Model  4165  of  the 
Series  58  bodies  with  our  potential  flow  and  boundary  layer  computer 
programs  without  any  iterations  to  add  boundary-layer  displacement 
thickness.  We  then  computed  the  drag  according  to  equations  (20)  and 
(21).  Comparison  with  the  measured  drag  reported  by  Gertler  (1950) 
is  given  in  table  2.  The  unreliable  predictions  obtained  with  equation 
(21)  do  not  necessarily  invalidate  the  equation.  They  may  be  caused  by 
inaccurate  values  used  for  the  boundary  layer  properties  at  the  end  of 
the  body.  This  region  is  the  most  inaccurate  for  the  boundary  layer 
calculation.  On  the  other  hand,  the  results  of  equation  (20)  are  not 
sensitive  to  the  inaccuracies  of  the  boundary  layer  calculation  in  this 
region  since  the  contribution  to  the  skin  friction  from  the  aft  end  i*s 
small. 


, 


As  indicated  in  figure  1,  the  parameterized  family  of  shapes  used 
for  Computer  Program  No.  1 consists  of  a nose  relation  and  a cubic  spline 
matched  with  continuous  slope  but  not  curvature.  Thus,  the  quantity  £ 
in  equation  (16)  was  set  to  zero,  the  quantity  a was  determined  so  as 
to  give  continuity  in  radius  at  the  match  point,  and  the  spline  was 
forced  to  have  the  slope  of  the  nose  shape  at  the  match  point  (see 
figure  5) . 

The  resulting  free  parameters  are  the  exponent  in  the  nose 

relation 


the  radii  at  the  spline  points,  except  for  the  tail  point  (this  radius  was 
fixed  at  zero) , and  the  slope  at  the  tail.  This  parametric  definition 
provides  a diverse  family  of  shapes,  including  some  with  blunt  noses 
(in  other  words,  those  with  a.  < 1 ) . 


Boundary  layer  transition  was  measured  to  be  about  seven  percent  of 
the  length  downstream  from  the  nose. 

trip  strip  five  percent  of  the  length  downstream  from  the  nose 
provided  a turbulent  boundary  layer. 


Starting  Shape,  Series  58-No.  4158, 
L/D  - 8.0  C - 0.0177 
°V 

Final  Shape, 
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3.2  Computer  Program  No.  2 

The  fluid  mechanical  calculations  for  Computer  Program  No.  2 are 
based  entirely  on  the  methods  used  by  Nakayama  and  Patel  (1974).  The 
potential  flow  is  calculated  by  the  method  of  Landweber  (1959),  and 
the  flow  field  is  constructed  with  doublets  distributed  on  the  axis  of 
symmetry.  The  governing  equation,  a Fredholm  integral  equation  of  the 
first  kind,  enforces  flow  tangency  at  the  body  surface.  It  is  solved  by 
iteration  with  numerical  evaluation  (by  Gauss  quadrature)  of  the  integral 
appearing  in  the  equation. 

The  turbulent  boundary  layer  calculation  method  in  Computer  Program 
No.  2 was  developed  by  Patel  (1974).  It  uses  power-law  velocity  profiles 
in  a momentum  integral  relation,  and,  for  auxiliary  relations,  it  uses 
an  appropriate  skin-friction  law  and  an  entrainment  equation. 

The  drag  is  calculated  in  Program  No.  2 according  to  the  method  of 
Nakayama  and  Patel  (1974),  and  the  explanation  given  in  the  reference  is 
repeated  here  for  convenience.  The  potential  flow  is  calculated  without 
boundary  layer  displacement  modeling  to  obtain  the  pressure  distribution 
for  the  boundary  layer  calculation.  After  each  step  in  the  boundary 
layer  calculation  over  the  last  ten  percent  of  the  body,  the  drag  formula 
of  Granville  (1953),  (a  modification  of  equation  (21)), 


C 


D 


7H+17 

8 


(22) 


is  applied  to  determine  a "local  drag  coefficient."  The  successive 
values  of  this  drag  coefficient  increase  up  to  some  point  and  decrease 
as  the  boundary  layer  calculation  approaches  the  tail  of  the  body. 
Instead  of  using  the  value  at  the  tail  end,  as  indicated  by  the  original 
Squire-Young  drag  method,  Nakayama  and  Patel  take  the  drag  coefficient 
simply  to  be  the  maximum  value.  Comparisons  of  the  drag  computed  in 
this  manner  with  experiment  are  given  by  the  authors  and  support  the 
acceptability  of  this  procedure  for  traditionally  shaped  bodies.  A 
recent  note  by  Patel  and  Guven  (1976)  reports  that  this  method  is  not 
reliable  for  bodies  with  thick  aft-end  boundary  layers  and  low  fineness 
ratios. 

For  Computer  Program  No.  2,  the  parameterization  of  the  shapes  is 
the  same  as  that  used  for  Computer  Program  No.  1,  except  that  continuity 
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through  the  second  derivative  is  enforced  at  the  match  point  between  the 
nose  relation  and  the  cubic  spline.  Thus,  the  full  nose  relation 

al 

r = ax  + Sx  (16) 

is  used,  and  the  quantities  a and  8 are  determined  by  the  continuity 
requirements . 
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4.  Calculated  Examples 

One  example  calculation  using  Computer  Program  No.  1 and  two 
calculations  using  Computer  Program  No.  2 are  presented  in  the  following. 
The  Reynolds  number  based  on  the  cube  root  of  the  volume  is  ten  million 
for  all  of  the  examples;  in  other  words,  Re^  = 10^  . Also,  the 
boundary-layer  transition  point  was  fixed  at  the  five  percent  body 
station  (X/L  = 0.05)  for  all  calculations. 

For  each  calculation  the  starting  and  final  shapes  are  illustrated, 
and,  for  the  final  shapes,  plots  of  the  boundary-layer  edge  velocity  'J  , 
the  planar  shape  factor  H , and  the  planar  momentum  thickness  Reynolds 
number  , as  calculated  by  the  hydrodynamics  part  of  the  computer 

programs,  are  presented. 

4.1  Calculation  No.  1 

For  the  first  calculation  we  used  Computer  Program  No.  1,  starting 
with  the  Series  58,  Model  No.  4158  body  shape.  The  spline  point 
abscissas  were  equally  spaced  along  the  x axis,  as  shown  in  figure  5. 
Therefore,  the  seven  parameters  varied  during  the  optimization  search 
were  the  exponent  of  the  nose  relation 

al 

r = ax  , (17) 

which  extended  from  x = 0 to  x = 1/6  , the  radii  at 
x = 1/6,  2/6,  3/6,  4/6,  5/6  , 
and  the  slope  at  the  tail. 

As  mentioned  earlier,  the  potential  flow  solution  involves  approxi- 
mating a surface  source  distribution  with  constant-strength  source 
panels.  For  Calculation  No.  1 we  obtained  the  paneling  by  equally 
dividing  the  axial  coordinate  from  the  nose  to  the  tail  and  by  equally 
dividing  the  polar  coordinate,  associated  with  a body  of  revolution, 
from  0 to  2tt  . For  both  directions,  the  number  of  divisions  was  chosen 
to  be  20. 
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The  calculated  drag  of  the  starting  shape  was 

C = 0.0177  , 

UV 

start 

and  the  measured  drag  was  reported  by  Gertler  (1950)  to  be 

C = 0.0180  . 

V 

The  resulting  minimum  drag  shape  is  shown  in  figure  5 together  with 
the  starting  shape.  The  calculated  drag  for  the  final  shape  was 

C = 0. 0144  . 

v 

final 

The  final  shape  is  characterized  by  a low  fineness  ratio  (3.8,  compared 
to  8.0  for  the  starting  shape)  and  a delayed  and  severe  convergence  to 
the  tail  point. 

Also  in  figure  5 we  compare  the  boundary-layer  properties  for 
the  final  shape  to  those  for  a Series  58  body  having  about  the  same 
fineness  ratio.  We  obtained  these  results  from  the  hydrodynamics 
part  of  Computer  Program  No.  1. 

We  note  here  two  features  of  the  velocity  distribution  shown  in 
figure  5.  The  dip  in  the  velocity  at  X/L  = 0.18  is  associated  with  a 
discontinuity  in  curvature,  which  occurs  where  the  nose  relation  is 
matched  to  the  cubic  spline.  For  the  calculations  with  Computer  Program 
No.  2,  for  which  curvature  is  continuous  everywhere,  irregularity  at 
this  match  point  did  not  appear.  The  second  irregularity  in  the  velocity 
distribution  in  figure  5 occurs  near  the  tail.  The  streanwise  panel 
length  used  for  the  numerical  approximation  to  the  surface  source  distri- 
bution for  the  potential  flow  calculation  is  larger  than  the  radius  over 
the  last  twenty-five  percent  of  the  body  length  (for  example,  at  X/L  = 
0.90  the  panel  length  is  five  times  the  radius).  We  believe  that  this 
coarse  paneling  is  responsible  for  the  velocity  irregularity  in  this 
region.  This  irregularity  did  not  occur  in  Calculation  No.  3,  where  the 
same  shape  was  analyzed  (as  the  starting  shape)  by  Computer  Program  No. 

2 using  a more  appropriate  discretization  for  the  numerical  solution  to 
the  potential  flow  problem  (see  section  4.3). 
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4.2  Calculation  No.  2 

For  Calculation  No.  2 we  repeated  Calculation  No.  1 with  the 
following  changes.  Computer  Program  No.  2 was  used,  and  the  spline 
point  abscissas  were  changed  to  decrease  their  spacing  at  the  front 
part  of  the  body,  as  shown  in  figure  6.  Thus,  the  seven  parameters 
varied  during  the  search  procedure  were  the  exponenent  of  the  nose 
relation 


1 

r = ax  + Sx  , (16) 

which  extended  from  x = 0 to  x = 0.075  , the  radii  at 
x = 0.075,  0.225,  0.425,  0.625,  0.824  , 
and  the  slope  at  the  tail. 

For  the  potential  flow  calculation  the  order  of  the  Gauss  quadrature 
formula  used  to  evaluate  the  integral  in  the  governing  equation  was 
chosen  to  be  24. 

For  this  calculation,  the  drags  for  the  starting  and  final  shapes  were 

C = 0.0180 

V 

start 

and  C = 0.0174  , 

v 

final 

and  the  starting  and  final  shapes  and  the  final  boundary-layer  properties 
are  shown  in  figure  6.  The  search  increased  the  maximum  thickness  of 
the  shape  somewhat  and  moved  the  point  of  maximum  thickness  forward. 

Table  3 summarizes  comparisons  of  the  first  two  calculations. 


5 

L 

Layer  Chat 
le  Same  Sta 
Program  N'o 
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4 . 3 Calculation  No.  3 

Since  the  results  of  Calculations  Nos.  1 and  2 are  so  markedly 
different,  we  thought  it  would  be  interesting  to  investigate  the 
importance  of  the  starting  shape  on  the  final  result.  For  this  reason, 
we  used  the  unorthodox  shape  obtained  as  a minimum  drag  body  by  Computer 
Program  No.  1 as  a starting  shape  for  Calculation  No.  3,  using  Computer 
Program  No.  2.  We  used  the  same  set  of  spline  point  abscissas  for  both 
calculations,  as  indicated  in  figures  5 and  7.  Also,  as  in  Calculation 
No.  2,  the  order  of  the  Gauss  quadrature  was  24  for  the  potential 
flow  calculation. 

The  results  of  Calculation  No.  3 are  shown  in  figure  7.  Note 
that  the  drag  for  the  starting  shape  is  considerably  higher  than  the 
value  computed  by  Computer  Program  No.  1: 

0.0144  for  Computer  Program  No.  1 
0.0186  for  Computer  Program  No.  2 


Also  of  interest  is  that  the  final  shape,  although  considerably 
thicker  than  that  for  Calculation  No.  2,  has  nearly  the  same  drag: 


0.0174  for  the  optimized  value  from  Calculation  No.  2 
0.0172  for  the  optimized  value  from  Calculation  No.  3 


10  Million 
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5.  Discussion  and  Conclusions 

The  optimization  calculation  using  Computer  Program  No.  1 produced 
a thick  body  (fineness  ratio  3.8)  with  a drag  19  percent  lower  than 
that  for  the  Series  58,  Model  4158  body,  a conventional  shape  with  a 
fineness  ratio  of  8 and  the  shape  used  to  initiate  the  optimization 
search.  According  to  the  check  cases  discussed  in  section  3.1, 
the  drag  model  used  in  Computer  Program  No.  1 underpredicts  the  drag 
for  low  fineness  ratios.  Figure  4,  based  on  the  information  in 
Table  1 of  that  section,  gives  the  discrepancy  between  calculated  and 
measured  drags  as  a function  of  the  fineness  ratio.  Although  the 
applicability  of  this  plot  to  other  shapes  is  not  established,  it  is 
useful  for  making  rough  corrections,  and,  if  applied  to  the  results 
of  the  optimized  shape  obtained  by  Computer  Program  No.  1 (see  figure  5), 
the  drag  is  increased  to 


C = 0.0175  , 

V 

which  is  nearly  the  same  as  for  the  starting  shape. 

We  expect  from  the  results  shown  in  table  1 of  section  3.1 
that  the  optimization  search  would  lead  to  a low  fineness  ratio.  It 
is  unfortunate  that  the  low  drag  calculated  is  apparently  the  result 
of  the  inadequacy  of  the  drag  model  used  and  not  the  result  of  a superior 
shape.  It  was  primarily  because  of  this  inadequacy  that  we  developed 
Computer  Program  No.  2.  We  note,  however,  that  the  optimization  method 
performed  satisfactorily;  it  successfully  found  a lower  drag  body. 

Computer  Program  No.  2 is  based  on  the  method  of  Nakayama  and 
Patel  (1974),  which  is  generally  regarded  as  the  best  state-of-the- 
art  drag  prediction  method  using  an  integral  boundary-layer  procedure 
(for  example,  see  Hess  and  James  (1975)). 

The  optimization  calculations  with  Computer  Program  No.  2 support 
the  conclusion  that  shaping  alone  cannot  significantly  reduce  the  drag 
below  that  of  the  conventional  shapes  exemplified  by  the  Series  58 
bodies.  In  fact,  the  two  optimized  shapes  calculated  had  drags  of 
only  five  to  six  percent  less  than  for  Model  4158  of  Series  58.  The 
conclusion  that  the  drag  is  insensitive  to  a shape  perturbation  from 
a streamlined  unseparated  body  when  the  boundary-layer  is  turbulent 
over  most  of  the  wetted  surface  agrees  with  previous  investigations; 
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for  example,  that  of  Hess  and  James  (1975).  We  note,  however,  that  for 
the  unusual  starting  shape  used  in  Calculation  No.  3,  the  drag  model 
predicts  a drag  significantly  higher  than  that  of  traditional  shapes, 
and  the  optimization  search  successfully  drives  the  shape  toward  a 
conventional  one. 

We  also  note  that  the  boundary  layer  shape  factor  H never 
exceeded  1.82  for  the  final  shapes  obtained  by  Computer  Program  No.  2. 

The  highest  value  occurred  at  the  tail  for  Calculation  No.  3 and 
corresponds  to  a value  of  2.59  for  the  planar  shape  factor  H , as 
indicated  in  figure  7.  Thus,  boundary  layer  separation  would  not 
be  expected  for  these  shapes,  and  the  inability  of  the  program  to 
predict  separation  and  calculate  the  associated  drag  is  therefore 
unimportant. 

If  we  had  used  more  accurate  methods  to  calculate  the  drag, 
we  may  have  been  able  to  show  more  significant  drag  reduction  by 
shaping  than  was  demonstrated  in  this  study.  The  example  calculations 
based  on  best  state-of-the-art  drag  prediction  method  using  an  integral 
boundary  layer  procedure  indicated  that  the  drags  of  bodies  with  all- 
turbulent,  unseparated  boundary  layers  do  not  differ  greatly.  However, 
the  drag  calculations  that  lead  to  this  conclusion  rely  on  the  formula 
of  Granville  (1953,  see  equation  (22)),  which  requires  an  evaluation  of 
the  boundary  layer  momentum  area  at  the  tail.  Since  standard  boundary 
layer  calculations  are  not  valid  near  the  tail,  it  is  difficult  to 
apply  this  formula.  Furthermore,  recent  observations  by  Patel  and 
Guven  (1976)  show  that  the  Granville  formula  gives  appreciable  error 
when  the  tail  boundary  layer  is  unseparated  but  thick.  The  breakdown 
of  the  boundary  layer  calculation  and  the  Granville  drag  formula 
results  from  neglect  of  radial  pressure  gradients.  In  addition,  thin 
boundary  layer  approximations  are  violated.  It  appears  that  more 
sophisticated  numerical  methods  (for  example,  finite-difference 
calculations)  must  be  used  to  calculate  the  drag  with  a high  degree 
of  accuracy.  These  methods  must  allow  for  the  following: 

(1)  The  boundary  layer  thickness  is  not  assumed  to  be  small 
compared  to  the  body  radius, 

(2)  The  pressure  gradient  normal  to  the  surface  is  not  neglected,  and 

(3)  The  viscous-inviscid  flow  interaction  is  adequately  represented. 
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We  believe  that  more  accurate  drag  calculations  incorporating  the  above 
items  should  be  developed.  Such  methods  would  be  useful  for  calculating 
the  turbulent  aft-end  region  of  laminar  flow  bodies,  as  well  as  all- 
turbulent  boundary  layer  bodies. 
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Appendix  A:  The  Potential  Flow  Boundary  Condition  Required  To 

Simulate  Boundary  Layer  Displacement 

There  are  two  methods  to  _ unt  for  the  displacement  of  the  boundary 
layer  when  obtaining  the  potential  flow  solution.  One  method  is  to 
displace  the  surface,  and  the  other  is  to  eject  fluid  from  the  surface 
(in  other  words,  to  prescribe  a surface  normal  component  for  the 
velocity  at  the  wall) . In  the  following  we  develop  an  equation  for  the 
latter  in  terms  of  the  boundary  layer  growth.  Important  quantities  are 
defined  in  figure  A-l.  Note  that  the  x and  y directions  refer  to  wall 
tangential  and  normal  direction,  respectively. 


Figure  A-l  Nomenclature 
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Two  velocity  fields  are  relevant.  The  first  is  that  the  velocity 
field  for  the  real  flow  (with  a boundary  layer)  for  which  the  x and  y 
components  of  velocity  are  denoted  by  u and  v , respectively.  For 
this  field,  the  velocity  vanishes  at  the  surface,  and  the  tangential 
velocity  at  the  edge  of  the  boundary  layer  is  denoted  by  U . The  second 
velocity  field  considered  is  a potential  flow  that  matches  the  real 
flow  outside  the  boundary  layer.  Its  components  are  denoted  by  u^  and 
v and  differ  from  u and  v in  the  boundary  layer  region  only. 

Note  also  in  figure  A-l  that  the  edge  of  the  boundary  layer,  a 
distance  6 above  the  surface,  is  not  a streamline.  We  will  consider 
some  streamline  just  outside  the  boundary  layer,  which  we  call  the  outer 
streamline.  Its  distance  above  the  surface  is  denoted  by  D . 

For  the  real  flow  field,  the  volume  flux  under  the  outer  streamline 


J u2iTr  dy  = J 2nUr  dy  - J 2tt(U  - u)r  dy  . 


Conservation  of  mass  requires  that 


so  that 


J 27TUr  dy  - — j 2tt(U  - u)r  dy  = 0 . 


For  the  potential  flow,  the  volume  flux  under  the  outer  streamline 


Qp  = J Up2:rr  dy  * j 2irUr  dy  - J 2tt(U  - Up)r  dy  • 


Conservation  of  mass  requires  that 


3— * 2irr  v , 
dx  op 

ro 

where  r = surface  radius 
o 

Vp  - normal  velocity  component  of  the  potential  flow  at  the 
0 surface. 


r 
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Thus , 


d_ 

dx 


J 2irUr  dy  - J 2tt(U  - u^)r  dy  = 2irrovp 


o o 

Subtracting  the  conservation  of  mass  equations  for  the  two  velocity 
fields  while  noting  that 


u = u for  y > 6 
P - 


gives 

6 6 
J 2tt(U  - u)r  dy  - J 2tt(U  - u^)r  dy 
o o 

We  assume  that  the  potential  flow  velocity  does  not  vary  significantly 
from  y = 0 to  y = <5  so 


d_ 

dx 


2iTr  v 

° P, 


u = 

[u  I 

i 

p 

l P) 

y = 6 


U , 


and,  therefore,  the  second  integral  on  the  left-hand  side  of  the  above 
equation  is  small  compared  to  the  first.  Thus,  the  normal  component  of 
velocity  at  the  surface  for  the  potential  flow  that  simulates  the  boundary 
layer  displacement  is  approximately 

5 

~ ~ 1 ^ f 27r(U  - u)r  dy  = -5—-  ~ (UA,)  , 


p 2?rr 
o o 


2nr  dx  'u“l' 
o 


where  the  boundary  layer  displacement  area  A^  is  defined  by 


^ = 2ir  J (1  - -jj)r  dy  . 


o 


